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Abstract 

In this work we propose an extended formulation for the interaction be- 
tween neutrinos and gravitational fields. It is based on the parametrized 
post-Newtonian aproach, and includes a violation of the universality of the 
gravitational interaction which is non diagonal in the weak flavor space. We 
find new effects that are not considered in the standard scenario for viola- 
tion of the equivalence principle. They are of the same order as the effects 
produced by the Newtonian potential, but they are highly directional depen- 
dent and could provide a very clean test of that violation. Phenomenological 

consequences are briefly discussed. 
PACS numbers: 04.80.Cc, 14.60.Pq 
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I. INTRODUCTION 



Despite the great success of general relativity to explain the gravitational interaction, 
this theory has proved to be very difficult to test in detail. There are several systems 
where gravitation is important, such as the very early universe, pulsars, quasars, black 
holes, and gravitational waves. However, gravitational fields in astrophysical systems can 
be considered as weak, even in the extreme cases of the neighborhood of a neutron star or 
at few Schwartzchild radii from a black hole. Usually the gravitational effects beyond the 
Newtonian level are very small and too tightly interwoven with other local physical effects 
to be clearly observed. 

The interest in these effects is twofold. They can shed new light on the character of the 
gravitational interaction, and they can provide invaluable information on some astrophysical 
systems, such as supernovas and neutron stars. With respect to the first point, two ques- 
tions arise naturally: the validity of general relativity as a description of the gravitational 
interaction, and the universality of this interaction. 

In fact, the strongest evidence of the universality of the gravitational interaction involves 
electrons, protons and neutrons, that is the members of the lightest family of matter fields 
in the standard model. Tests of the weak equivalence principle for these particles include 
laboratory experiments of the Eotvos-type, which measure the gravitational acceleration 
of macroscopic bodies. They state that gravity accelerates all macroscopic objects at the 
same rate to an accuracy of one part in 10^^ IQ. The experimental limits for the universality 
between matter fields and gauge fields are weaker. For example, the supernova SN1987A gave 
the opportunity of a direct comparison between the transit time for photons and neutrinos 
traversing the same path in a gravitational potential 0(r), which leads to limits on the 
violation of the weak equivalence principle by massless particles of the order of [7^ — 7i/| ^ 
10^'^, where 7 is the PPN (parametrized post-Newtonian formalism) parameter for the scalar 
potential 0. 

It is more difficult to obtain observational evidence on the gravitational coupling of the 
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heavier families, except for the case of the kaon system where the bound I0A7I < 2 x 10~^^ 
0] has been set. Information of this kind can also be acquired from the propagation of 
neutrinos in a medium with a gravitational background. As is well known, ordinary matter 
affects the neutrino propagation in a flavor dependent way and, under favourable conditions, 
large transformations of one neutrino flavor into another can take place, even for small 
mixing between the mass eigenstates. The implications of this mechanism in astrophysics 
and cosmology has been extensively examined during the last years. A similar resonant 
enhancement could be induced by a non-universality in the gravitational interaction of the 
neutrinos, even if they are massless. Violations of the equivalence principle of the order 
of 10"^'' cannot be ruled out in this context, and in fact the observed deficit in the solar 
and atmospheric neutrino fluxes have already been interpreted as a positive signal of this 
violation An interpretation in terms of a violation of the Lorentz invariance is also 

possible 0, but this could be included within the violation of the equivalence principle 
scheme in the case of constant fields 0]. Note that neutrinos are unequaled as test particles 
for probing the gravitational field. Because of the smallness of their interactions, the level 
of accuracy that can be achieved with them is several orders of magnitude better than in 
any other previous test. 

The present work is partially motivated by the above considerations and develops a gen- 
eral framework for analyzing the possible flavor dependence of the gravitational interaction. 
Our approach is a generalization of the one proposed in Refs. @J^, where only the effects 
due to the scalar gravitational potential was considered. In this way, we find new contri- 
butions to the oscillations that are of the same order of magnitude than the terms involving 
the Newtonian potential, and we extend the analysis to the next PPN order, which includes 
contributions generated by the angular momentum of the gravitational source. In contrast 
with the scalar potential contribution, this gives place to new highly anisotropic effects, 
which in principle could be verified in several astrophysical systems. These new effects pro- 
vide a more precise and characteristic signature for a possible violation of the equivalence 
principle. 
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Whereas the integer spin fields can be consistently described in a curved space-time, 
the half-integer spin fields need to be defined with reference to a locally inertial frame at 
each point of the space- time [Q. If we are considering several neutrino fiavors, we have 
to extend the usual construction by defining each gravitational fiavor in its own inertial 
frame. Furthermore, we can introduce a possible violation of the equivalence principle and 
assume that these frames are not necessarily related by Lorentz transformations, and thus 
could be physically non-equivalent. At each of the orthonormal frames the gravitational field 
is supposed to have the structure given by the PPN formalism, which provides a general 
account of the possible deviations from the Einstein theory. 

These assumptions lead us to a generalization of the standard scenario for the violation of 
the equivalence principle (VEP). With our present theoretical understanding it seems that if 
we want to keep the spinor structure, the theory does not correspond to a metric one which 
can be settled on a consistent basis. More precisely, our approach should be considered as 
a phenomenological one, which allows us to search for possible signatures of fiavor-changing 
effects associated with a violation of the equivalence principle. In the following we do not 
consider the dynamics of the gravitational fields but we still use a manifold as the space-time 
framework, i.e., although the metric is not defined we assume that the space-time is well 
defined. 

II. NEUTRINOS IN A GRAVITATIONAL BACKGROUND 

Given the frame the equation for a (massless) neutrino is the Dirac equation in a 
curved space time: 

7"v;^(9^ + ^r^)^ = o, (1) 

where the connection is 

r, = -iKfv,r^V/3, (2) 
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with (Jq,/3 = f [7a, 7/3]- All covariant derivatives are metric, so we are neglecting small torsion 
effects . Other possibilities for the equation involve generical couplings with the curvature 
T0| , but these terms are highly suppressed in the usual astrophysical situations by the small 



values of the gradients of the gravitational fields as compared to the neutrino momentum. 



Explicitly, the linearized Dirac equation for a static gravitational field reads |]TT 



^I'd, - {/loo, 7 ■ V} - ^ 7^9^} + {/loi, 7°5^} + ]p\'=^d^h'^s^^ vl>^ = , (3) 

where the h^'^ fields are defined by g'^^ = rj^^ + h^^ and = ^e^^^a^K The spatial derivatives 
of the gravitational potentials are proportional to the inverse of their characteristic variation 
length, dih^^ oc . In astrophysical systems the spatial dimensions and the energy of 
the neutrinos render L ^ X^^ where oc p^^ is the neutrino wavelength. Accordingly, 
we can neglect the terms with spatial derivatives of the gravitational fields, including the 
spin contributions. This approximation can be justified on a more general basis by using 
a geometric optic-like expansion of the Dirac equation in powers of the parameter X/L 
||12|| . The main effects of the neglected terms, which include a chirality transition induced 
by gravity, are independent of the equivalence principle violation, and have been already 
analyzed in Ref. []rT|. In the above approximation Eq. (|^) reduces to 



iao*. = H^, , (4) 

with the Hamiltonian given by 

H = -i7o7i[(l - ^/i°°)(9i - ^hijdj] - ihoidi , (5) 
where all metric dependent terms are assumed to be slowly varying functions of the position. 



In an astrophysical scenario we can work within the framework of the PPN theories [jT3 
The assumptions for constructing the metric in the PPN formalism involve virialized sources 
such that ^ ~ w^, where the quantities M,R, and w represent estimations of the order of 
magnitude of the mass, distance and characteristic (average) velocity of the source. 
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The metric is the Minkowskian one plus source dependent perturbations. The latter 
have the correct tensorial character and dimensions, falling at least like at infinity. 
This metric is in general given by 

Ko = 2iU + O(w^) , (6) 
hoi = —^iVi - ^A^Wi + («2 - ^«i)vif/ - a^yjU.i + C>(w^) , (7) 
= 2^U6^j + TUij + O(w^) , (8) 

where the potentials are 

U=[P^^, (9) 
J I r — r' I 

_ f p{r'){ri - r'i){rj - r'j) rfV 

~ J I r - r' |3 ' ^ ^ 

J I r — r' I ' 

/■ p(r')(w(r') ■ (r — r'))(r,- — r') (i^r' ^ ^ 

= J \r-J\^ • (^2) 

Here p(r) is the density of mass and w(r) is the velocity of the source of the gravitational 
field. We are using a system of unities where G = h = c = 1. We keep each term of the 
same order in the PPN expansion because we are interested in ultrarelativistic neutrinos. 

In the particular case of a very confined and distant source, the expressions in Eqs. 
(D-(IO) can be approximated as follows 

V,^w,U+0 (i,) , W, « + O (-1) , (14) 

where Xj are the components of R. 

Up to order w'^ the adimensional parameters of the expansion are 7, 7', Ai, A2, F, v , 
«! and a2- In Einstein gravity we have ai = 0:2 = T = , 7 = 7' = Ai = A2 = 1 and v is 
irrelevant. The parameters ai and 0:2 are null if the theory is Lorentz covariant, but if there 
is a preferred reference frame, characterized by a velocity v, they should be non null. In 
general the characteristic frame velocity can be a gravitational fiavor dependent quantity. 
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The parameter ai can be fixed to be 7Ai + A2 — 47 — 47', while 0:2 is an independent 
parameter up to this order. 



where E is the energy eigenvalue and we keep only first order terms in the metric /i^jy. 

III. NEUTRINO OSCILLATIONS FROM A NON-UNIVERSAL 
GRAVITATIONAL INTERACTION 

The purpose of our work is to state a meaningful phenomenological basis to discuss a 
possible violation of the equivalence principle. The breakdown of the universality of the 
gravitational interaction raises the question of the complete consistency of the underlying 
theory, but we will put aside this issue. Here, we are concerned only with a model that 
describes the different phase shifts associated to a possible flavor dependence of the gravita- 
tional interaction, its consequences, and the possibilities of its detection in an astrophysical 
framework. We assume that the metric associated with each family are very close to each 
other and to the Einstein theory predictions. The PPN coefficients depend on the flavor and 
are assumed to be diagonal matrices in the gravitational flavor basis. The more general case 
where these coefficients are arbitrary matrices would imply that there exist no local inertial 
frames where the neutrinos can be deflned, and therefore their spinorial nature would not 
be clear. We do not consider this situation. 

Following this approach, each of the n neutrino flavors is deflned in a different orthonor- 
mal frame {V^"^}, a = l,...,n, which form non-equivalent bases for the tangent space. 
Accordingly, we write the metrics as 



The dispersion relation implied by the Hamiltonian in Eq.(^) is 




(15) 



a 



(16) 



From 




neutrino gravitational flavor can be 



approximated by 
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E^=p{l + h% p. - U,,p.p, - (7'" + 7")f/) ■ (17) 

In the PPN approximation the coordinate system is generally fixed to give F = and 7' = 1, 
(where the last equality is equivalent to the definition of the Newton constant,) but in our 
case such a coordinate fixing could be done only for one metric at the expense of the others. 
Therefore, we are setting these parameters at the usual values for the first gravitational 
family, while the others are left as free parameters. This is an important point as we will 
see in what follows. 

If the parameters are family dependent, then distinct neutrinos will undergo different 
phase shifts when passing through the same sector of the space. The phase shift differences 
become observable when the particle basis that diagonalize the weak and the gravitational 
interaction are not the same. In this context, the effects of a universality violation at the 
level of the scalar potential U has been already considered in Ref . 0] , where it is shown that 
such a violation leads to neutrino oscillations. The factor (1 + 7*^) of their model should be 
replaced by (7"^ + 7"^) in our case. We extend this model by including the effects of the Uij 
potentials that are of the same order of magnitude as U, together with the PPN structure 
expanded up to w^, with the corresponding generalized universality violations. In this way, 
we are taking into account not only the coupling of the neutrinos to the mass of the source of 
the gravitational field, but also the coupling to its quadrupolar distribution and its angular 
momentum. As we will show, these last interactions produce highly directional and very 
characteristic effects. 

In what follows we examine the neutrino oscillations induced by a non-universal grav- 
itational coupling. To find out the main features of this phenomena we will consider two 
neutrino fiavors. We assume that the gravitational fiavor basis is related to the electroweak 
basis through a unitary transformation U, characterized by a mixing angle 6g: 











cosOg -sin 6g ^ 












^sin 6g cosOg ^ 





The equation for the neutrino evolution is 
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(19) 



where ?i is a diagonal matrix in the gravitational flavor basis, whose eigenvalues are given 
by Eq. (IT^). After discarding an irrelevant overall phase, we have 



(20) 



with 



d 




_ Ao 


-cos26g sin 26g ^ 


(•^ \ 


I— 
dt 




2 


^ sin 26g cos26g ^ 




E2 - 




= E [6hoiPi - 6T UijPiPj - 





(21) 



where E = p is the neutrino beam energy, and 



Shoi — - /iQj 



7 1 1 

--SAiVi - -SA2Wi + {Sa2 - -Sai)viU - Sa2VjUji . 



(22) 
(23) 



For a constant gravitational field the survival probability is 



1 - sin^ (2^3) sin^ — 
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(24) 



where L = t — \s the distance travelled by the neutrino from the production point. 
This clearly shows that oscillations will appear whenever there exists a non null mixing 
angle induced by flavor dependent gravitational interactions. These oscillations have a 
characteristic length given by 

27r 



I An 



(25) 



In contrast with vacuum oscillations induced by a mass difference, where = |^ is 
proportional to the energy, the effect we are considering has oscillation lengths proportional 
to which makes this phenomena suitable to be observed in the case of high energy 

neutrinos. Note that even though the overall sign of the gravitational potential is irrelevant 



for the oscillations, the relative signs among the parameter differences are very significative. 
The Eq. (|2^ leads to the following averaged survival probability 

< P(z/e Ue) >= ^ (l + COS^^Og) ■ (26) 

As is well known, neutrino oscillations in matter are qualitatively different from the 
oscillations in the vacuum. This is because the interaction of the neutrinos with matter 
modify their dispersion relations. Neutral current interactions are flavor diagonal and can 
be ignored, as long as we do not consider sterile neutrinos and neutrinos are not part of 
the medium. In general, this will not be true for the charged current interactions. The 
forward scattering amplitude is not flavor diagonal in this case, and depends on the leptonic 
content of the matter. The above gives place to important consequences such as the MSW 
effect. In general, we expect that a non universal gravity will also affect the electroweak 
Lagrangian, introducing a number of unknown coefficients, but the combined effect should 
be of the order U Gp and therefore is highly suppressed. 

If electrons are the only leptons that are present, then the matrix 

(i o \ 



bjt) 



(27) 



has to be added to the second term of Eq. (120). Here, be{t) = ^/2GFNe{t), with Ne{t) 
denoting the electron density. The resulting Hamiltonian 7i(t) can be diagonalized at every 
moment by introducing the instantaneous flavor basis, defined in an analogous way to Eq. 
(|1|), with Og ^ Omit) and 

sin2^„(t) = ^° . (28) 

^(Ao cos29g - he{t)Y + (Ao sin2^g)2 



In the adiabatic approximation, the average survival probability is given by the following 
formula 

< P(z/e ^ Ve) >= ^ (1 + cos2^g cos20„(to)) • (29) 
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which reduces to Eq. ( p6D when be(t) = 0, so that 9m(t) = 6g. The use of the adiabatic 
approximation is justified whenever 



R) 



dN,{t) 



dt 



I A 

° cos26'o 



(30) 



There also exists the possibihty of a resonant conversion when the diagonal elements of 
the full Hamiltonian vanish, i.e. when 



V2GFNe{tR) = Ao cos . 



(31) 



This mechanism can totally change the fiavor, independently of the value of the mixing angle 
9g, but its efficiency depends on the adiabaticity of the process. 



IV. PHENOMENOLOGICAL EFFECTS 

There is an inherent uncertainty in the potentials in Eqs. (^-(10), because arbitrary 
constants can always be added without changing the physics, as far as the effects associated 
with a violation of the equivalence principle are not involved. In Einstein theory these 
constants can be eliminated by a coordinate transformation. Instead, in the PPN expansion 
the coordinate system is fixed, so possible uncertainties arise from the very distant unknown 
mass distributions. We will restrict ourselves to the most important known near sources for 
the potentials, leaving aside the problem of the very distant ones. In any case, the distant 
sources would only produce significant isotropical effects, and thus could only affect the 
definition of U and the diagonal part of Uij (which in the formula for Xg can be absorbed 
into U). 



A. Solar neutrinos 

The sun is subject to a gravitational field that has several sources. The main sources are: 
our galaxy, the Virgo Cluster, and the Great Attractor. The more important contribution to 
the potential U comes from the Great Attractor gravitational field, with small perturbations 
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due to galactic clusters, and the galactic and solar fields. In consequence, this potential can 
be approximated by a constant of the order of 10~^ HHHl- The effect of this potential 
regarding a possible universality violation has been already analyzed in Ref. Resonant 
flavor oscillations are consistent with the observed deflcit of solar and atmospheric neutrinos 
for [7(5(7 + 7') — 10^^^. The effect is isotropic and to account for the more recent data 



requires a three-neutrino mixing scheme |15 



In the PPN approximation the contributions of the Uji potential are in principle consid- 
ered of the same magnitude as U. In the case of a distant source in the z direction, we have 
Uzz ~ U. However, the components Uxx, Uyy and U^y are proportional to {A6)'^U, where 

is the angular size of the source, while U^z and Uyz are of the order of A6U (in fact, 
they are proportional to the center of mass distance to the z axis). Considering that the 
Great Attractor is a rather extended object with an angular size of the order of 10^^ ||16|| , 
we see that in the case of the sun there are only three relevant types of Uji contributions: 
those coming from our galaxy, which are of the order of 10~^, a longitudinal component 
from the Great Attractor, of order of Uzz — f/ — 10^^, and transverse-longitudinal compo- 
nents also produced by the Great Attractor, of the same order as the galactic contributions, 
Uxz — Uyz — 10^^. The contributions of Vi and Wi are roughly proportional to source 
velocity w times U, which is of the order of (10"'^ — 10^^) f/. Their directional effect has a 
dipolar structure and can be assumed to be one or two orders of magnitude smaller than 
the dominant ones. 

Therefore, a violation of the equivalence principle would be characterized by three main 
effects manifested as flavor oscillations: an isotropic effect {U ^ 10~^), and two additional 
anisotropic effects {Uzz — 10^^, Uji ~ 10^^). If we assume that the differences of the PPN 
parameters due to the flavor dependence are all of the same order, the most signiflcant 
directional effect is given by a quadrupolar contribution due to Uzz- This effect could be 
of the order of the dipolar one originated by the elliptical orbit of the Earth, but the 
latter only depends on the eccentricity of the orbit (perigee: RA=18:48, DEC=-23:27 in 
equatorial coordinates), whereas the gravitational one depends on the energy of the neutrinos 
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and their direction with respect to the Great Attractor. The approximate position of the 
Great Attractor center in galactic coordinates is / = 325°, b = —7°, v = 4882 km s~^, 
or in equatorial coordinates is RA=a=16:10, DEC=(i=-60°10'. In ecliptic coordinates the 
aphelion position is RA=18:48 (282°), DEC=0°, whereas the Great Attractor center is at 
RA=A= 16:52 (253°), DEC=D=-38°17'. Both axis differ in approximately 30° in RA, and 
therefore the effects could be discriminated. Taking into account the dominant contributions, 
due to U and Uzz, and assuming U ~ Uzz, the neutrino wavelength in vacuum becomes: 

X = ^ (oo\ 

' EU\6Tcos^Dcos^a-A) + {S'y' + S'y)\' ^ ' 

where a is the right ascension of the sun in ecliptic coordinates at a given time. 

An easily visible consequence of the gravitational contribution is a breaking of the reflec- 
tion symmetry of the neutrino flux with respect to the aphelion-perihelion axis of the Earth 
orbit. This symmetry is characteristic of the scenarios which do not consider the violation 
of the universality of the gravitational interaction, with the only exception given by a pos- 



sible interaction between the solar magnetic field and the neutrino magnetic moment |]T7 



In general, the scenarios usually considered would yield different neutrino fluxes for Earth 
positions separated by six months. Otherwise, the gravitational contribution has the same 
sign for these two positions. 



B. Atmospheric neutrinos 

The dominant contributions in this case are the same as those already considered in the 
previous section. The main difference is originated by the Earth rotation which gives place 
to diurnal neutrino flux variation. This situation can be described more appropriately by 
means of azimuthal coordinates In terms of these coordinates the neutrino direction 
can be written: 

= (sin 9i, cos z^, sin 9^ sin Zj,, cos 9,^) , (33) 
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where is the zenithal distance and Zy is the azimuthal angle of the incident neutrinos. 
Similarly, for the Great Attractor position we have: 

dcA = (cos V? sin (i — sin cos d cos r, cos d sin r, sin ci sin (/? + cos d cos (p cos r) , (34) 

with being the Great Attractor declination angle and if the observatory latitude. The 
parameter r is a — t^, where a is the right ascension of the Great Attractor and tg is the 
sidereal time. According to Eqs. (^TD and (^51) the oscillation wavelength depends on 



{pu-dGA^ ■ If there were no violation in the universality of the gravitational interaction, 
then the neutrino flux would only depend on the zenithal distance 6,y. Here we have an 
additional dependence on the azimuthal angle z^, and the time r, which implies a diurnal 
variation of the flux. For instance, if we focus our attention on zenithal neutrinos {61, = 0), 
for the Kamiokande site {(f = 36.5°) we have: 

2n 



A 



q 

EU 



.5 + .4 cos 6T + (^7' + ^7) 

This is the effect at a given time r on the vacuum wavelength. It can be traced on with the 
dependence of the total flux integrated over r and z^. 



(35) 



C. Neutron star kicks and rotation 

Resonant neutrino oscillations influenced by the magnetic field in the early stage of a 
neutron star have been proposed as a possible mechanism to explain the observed proper 



motion of pulsars [|l^]. But to do this excessively strong magnetic fields are required ||20|| . 
Furthermore, the origin of the high angular velocity of the pulsars is also not completely clear 
||2T| . In the present discussion, the vectorial potential h^i manifests itself as an anisotropic 
perturbation in the shape of the resonant surface that is located between the neutrinospheres 
of the different flavors. This anisotropy could be the source of both the angular and the 
traslational accelerations, and thus could simultaneously give an explanation for the observed 
spins and proper motions. 
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The potential hoi contains two main contributions, both produced by the star. One of 
them, given by {6a2 — ^6ai)viU — 6a2VjUji , is relevant when a preferred frame exists, and 
therefore the coefficients a are non null. Their action is analogous to the one produced by 



a strong magnetic field [0,^], which can generate a translational kick in the movement of 
the star during neutrino emission. Taking reasonable values for the potential and a frame 
velocity of the order 10~^, the resulting gravitational effect has the correct magnitude to 
explain the observed kicks. The other contribution to hoi originates from Vi and is caused by 
the star rotation. It induces an angular acceleration during neutrino emission because the 
resonant surface depends on the neutrino angular momentum. The sign of this acceleration 
is determined by the relative sign between (^7' + ^7) and 6A1. 

Since the gravitational fields in the interior of a neutron star are relatively high, U ^ 1, 
the PPN approximation is not very good in this case. A more complete treatment should 
include higher order terms in the expansion, and hence more parameters. 



V. CONCLUSIONS 

In this work we have formulated a generalized scenario for neutrinos propagating in the 
presence of gravity. To incorporate the violation of the equivalence principle each neutrino 
fiavor has been defined in a different orthonormal frame at each point of space-time. This 
is a natural prescription if we are giving up the metric but still want to retain the manifold 
structure. The violation is parametrized by a generalized PPN expansion. 

In this way we have developed a VEP scenario for neutrino oscillations, which leads to 
several new effects. The most relevant is the one due to the potential Uji, which can be of 
the same order of magnitude as the Newtonian contribution. In the PPN metric theories 
this potential is irrelevant because it can be set to zero by a coordinate fixing. This is not 
the case in the present approach, where there is more than one metric. In addition, other 
potentials could be relevant for the early stages of neutron stars, giving place to angular or 
translational accelerations. 
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An interesting feature of the effects we have discussed is that they have a very character- 
istic directional dependence. These effects should be present in the solar and atmospheric 
neutrino fluxes if their anomalies are related to a VEP scenario. In such a way, these phe- 
nomena could provide relevant information on this violation. Accordingly, the solar neutrino 
flux must change along the orbit of the Earth in a way that clearly differs in the angular 
and energy dependence from the geometrical effect due to the orbit eccentricity and its 
consequences on the ordinary vacuum oscillations. In contrast to the mass-mixing scenario, 
the seasonal effect we consider could also appear in the case of resonant transitions in the 
sun, due to the variation of the position of the resonance region. In the case of atmospheric 
neutrinos, a detailed analysis of the daily and zenithal angle dependence of the flux, could 
reveal the effects we are interested in. Further analysis along these lines are in progress. 
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